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Abstract 

Experimental studies of meson production through two-photon fusion in inelastic electron- 
nucleus scattering is now under way. A high-energy photon radiated by the incident electron 
is fused with a soft photon radiated by the nucleus. The process takes place in the small-angle- 
Coulomb region of nuclear scattering. We expound the theory for this production process as well 
as its interference with coherent-radiative-meson production. In particular, we investigate the 
distortion of the electron wave function due to multiple-Coulomb scattering. 
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I. INTRODUCTION 

The PrimEx Collaboration [lj aims at measuring electromagnetic properties of pseu- 
doscalar mesons through the Primakoff effect, in 11 GeV/c electron-nucleus-inelastic scat- 
tering. Experiments of this type were suggested some years ago by Hadjimichael and Fallieros 
2!, and they also elaborated their theoretical description, within the Born approximation. 
Here, we shall develop, in a Glauber approach [3], a full-fledged theory for this process. We 
are in particular concentrating on the effects of Coulomb scattering of the electron and on 
the interference of the two-photon-fusion amplitude with the coherent-pion-photoproduction 
amplitude. The Primakoff effect on proton targets at JLab energies has been studied by 
Laget [4]. 

The kinematics of the electron-nucleus-pion-production reaction is defined through 

e-(h) + A(pO e'(k 2 ) + 7T°(K) + A(p 2 ). (1) 

Our analysis is carried out for high energies and small transverse momenta, meaning small 
compared with the longitudinal momenta. In addition, the momentum transfers to pion 
and nucleus must be in the Coulomb region, leading to further severe restrictions. In the 
PrimEx experiment typical energies are; for the incident electron E(kx) = 11 GeV, and for 
the scattered electron E(k 2 ) = 300 MeV. This implies an energy of 00 (K) = 10.7 GeV for 
the pion, as for the virtual photon initiating the pion-photoproduction process. The PrimEx 
Collaboration also plans to investigate production of 77 and 77' mesons. Our formulae are 
valid also in those cases but we shall present numerical results only for pions. 

The electron-nucleus production amplitude is a sum of two amplitudes; the two-photon- 
fusion amplitude M 2l , and the coherent-photoproduction amplitude A4 7 . The normaliza- 
tions are chosen so that the cross-section distribution takes the form 



da _ 1 



\M 2l + M 1 \\ (2) 



d 2 /W 2 A ± dfc 2 || 32(27r) 5 k 1 E(k 2 )K ll M 2 A 

with K\i the component of pion momentum along the incident ki direction. We neglect the 
recoil energy of the nucleus in comparison with its mass. 

The structure of the cross-section distribution at small-transverse-momentum transfers is 
essentially determined by the photon-exchange propagators. There are two such propagators; 
one in the variable k 2 ± and one in p 2± . Each of them exhibits a Primakoff-peak structure. 
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II. POINT-LIKE-NUCLEUS TARGET 



We start with the two-photon-fusion matrix element. In this initial calculation, the 
nucleus is treated as a point-like particle of zero spin. Moreover, we ignore the Coulomb 
phase of the electron wave function, which is generated by Coulomb multiple scattering of 
the electron. In other words, only the Born, or tree, diagram is considered. It is displayed 
in Fig. [U 

fci fc 2 



K 




FIG. 1: Graph describing pion production in the Coulomb field of a nucleus in inelastic-electron 
coherent-nucleus scattering. 

Naming the particle momenta as in Fig. [TJ the overall four-momentum conservation reads 

h + pi = k 2 + p 2 + K. (3) 
The four-momenta of the intermediate photons are defined according to 

fci = k 2 + q, (4) 

Pi = P2 + Q, (5) 
K = q + Q. (6) 

The vertex functions of the photon-fusion diagram are 

M{h -> k 2 + q) = M(fc 2 )ie7 M M(fci) (7) 
M( Pl ^p 2 + Q) = -iZe( Pl +p 2 y, (8) 
M{q + Q^K) = — g^e^ pu q p Q a . (9) 

TTl n 

The coupling constant g nll of Eq.([9]), with eoi23 = 1, is fixed by the anomaly of the axial- 
vector-current divergence. Taking its numerical value from the pion-decay rate 

r(7r° -> 77) = \ixa 2 gl n m^ (10) 



gives g irri = 0.0375. 

The complete matrix element for the photon-fusion diagram of Fig. [1] becomes 

it 2 1 1 

M 2l = ie(-iZe)( — g^){-i) 2 — • u(k 2 )^u(k l )^ vpa q p Q (7 (p l + p 2 f . (11) 
Tn^ q 

This expression can be simplified. Referring to Eq.flH]) we may write 

Pl +p 2 = 2 Pl -Q, (12) 

where the Q-term contribution to Eq. lflTj) vanishes due to gauge invariance, manifested by 
the anti-symmetry property of the e factor. Evaluating the matrix element in the lab. system, 
where p\ = (Ma, 0), we get 

M 2l = iZ—g^lMA \ ■ u(k 2 )j ■ (q x Q)u(/ci). (13) 
m v q z Q z 

So far, no approximations have been made. In the following reduction one must be careful 
since the vector q = (q_i_, q\\) has a longitudinal component q\\ = 10.7 GeV/c, and large, but 
a transverse component q± with a length rather measured in MeV/c. 

Expressed in two-component form the current-matrix element of the three- vector 7 = (7 fc ) 
can be approximated as 



u 



(*2)r*(*i) = (£i£ 2 ) 1/2 d 



1 i 
2ki + ^-k 2± - — (k 2 _L x a) 

£>2 &2 



Ci- (14) 



Here, we have used that the electron energies are much larger than the electron mass, 
Ei,E 2 ^> m e , and that for the scattered electron the transverse momentum is much smaller 
than the longitudinal one, k 2 ± -C k 2 n. 

Now, if we take advantage of the fact that all transverse momenta are small, q±, Q± <C q\\, 
that Q11 <C q», and are ready to neglect small terms of order E 2 /Ei, then the two-fusion 
matrix element simplifies to 

g4 E 1 1 [ 

M 2l = iZ—g^ 1 2M A \/l^E 2 -i- • — • C] (q± X Q ± ) -k x - q ± • Q ± icr k x Ci- (15) 
m n hj 2 q 1 J 

We observe that the spin- independent term is maximal when qj_ is orthogonal to Qj_, whereas 

the spin-dependent term is maximal when q^ is parallel to Qj_. 

Next, we put e 2 = Ana in the matrix element of Eq. (115j) . square it and take the proper 

average over the electron-spin states. Inserting the result into Eq.©, gives the cross-section 

distribution 

da 1 \Za 2 E 1 q ± Q±l 2 



with the double- Primakoff structure clearly exhibited. Observe, that in the denominators, 
q 2 and Q 2 are still squares of four- vectors. It will susequently be shown that in the kinematic 
configuration studied here this expression simplifies to 

dor i 



d 2 k 2± d 2 K ± dk 2ll nKu 



Za 2 q ± 
m n q]_ + m 



Q- 



Qi + Q 



(17) 



with Q\\ = m 2 j2K\\. 

In the following we shall, when no confusion is possible, use the notation .M27 also for 
the 2x2 production matrix, of which Eq.(TT5l) is the matrix element. 
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III. IMPACT-PARAMETER REPRESENTATION 



We shall eventually need to introduce the Coulomb phase of the high-energy electron. 
This will be done in the eikonal approximation. For this task we need the coordinate-space 
description of the two-photon-fusion amplitude. In photoproduction the incident photon 
combines with a virtual photon radiated by the electric field of the nucleus to form the 
produced pion 5j. In the present case there are two virtual photons involved; one radiated 
by the nucleus, the other radiated by the electron. Hence, two electric fields, or potentials, 
are needed. We start with the nucleus-radiation process. 

The kinematics is defined in Fig. 1. The values of the longitudinal momentum transfers 
are fixed and it is important to understand their sizes. In the energy-conservation constraint, 

^(kj) + M A = E 1 (k 2 ) + w(K) + E A (Q), (18) 

we expand in powers of the longitudinal momenta. Also for the scattered electron the lon- 
gitudinal momentum is much larger than the electron mass and the transverse momentum, 
so the expansion is possible also there. For the nucleus, however, the recoil kinetic energy 
is so tiny it can be neglected altogether. We obtain 

Q, - ^. (19) 

Since Kn = q\\ +Q\\ and Q\\ extremely small we may replace Kn by q\\ if we like. Numerically, 
at Kn = 10.7 GeV/c, the longitudinal momentum transfers for 7r°- and 77-meson production 
become 

Q^m^o) = 0.85 MeV/c, (20) 
Q\\{m v ) = 14.0 MeV/c. (21) 

The value for eta production is large on a nuclear scale. 

The photon radiated by the nucleus has four-momentum Q and we may put 

-Q 2 = Qi + Q|, (22) 

since Qo is negligible in comparison with the space components of Q. The action of the 
nuclear photon can be described by a vector field, V^r), which following Eq. ffT5]) has the 
structure 

Fa(q) = qT+q 1 = / d3re ~ iQ ' r VA(r) ' (23) 
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with Q = K — q. The index A indicates that the field is associated with the nucleus. 
The potential V^r) is exactly the one that appears in real photoproduction For a 
point-like-nucleus-charge distribution its analytic form is 

-1 s 



Va(t) 



with r = (s, z). 

Integration along the z-direction produces a vector-profile function, 



(24) 



with 



Fa(Q) = ^ J d 2 b e~^ h r A (b,Q||), 



r A (b,Q||) = - [bQuKxibQu 



(25) 



(26) 



and q, with q « K», the momentum of the photon radiated by the electron. The function 
Ki(z) is the modified Bessel function. 

The radiation by the high-energy electron is similarly described by a vector field, as can 
be seen by rewriting the propagator 



q 



k 2 \\ 



m 



(27) 



This structure leads to an energy- dependent electron potential. We pull out the energy- 
dependent factor —ki/k2z and define 



F e (q) 



m 



with q = ki — k2. The index e indicates that this vector field refers to the high-energy 
electron. Inverting Eq. (|28p . we get an expression for the corresponding potential. In analogy 
with that of Eq. f )24|) . and with r = (s, z), the electron potential reads 

1 s 



d 3 r e iq ' r VJr) 



(28) 



V e (r) 



Airi 



s 2 + 



k 2 



mi 



3/2 ' 



(29) 



The ^-dependence enters multiplied with an enormous scale factor. 
Defining the profile function through 



F e (q) 



2tt 



d 2 6e^' b r e (b,g,|), 



(30) 



and integrating over the z- variable in Eq. (128]) yields 

r e (b, 9 ,|) = ^ \ym e bKi{m e b)^ . (31) 

Due to the scaling of the z- variable, the longitudinal-momentum transfer enters in the com- 
bination m e q\\/ki ~ m e . 

We conclude from expressions ( 1261) and ( |3T|) for the profile functions that the cut-off 
in impact parameter b is 2Kn/m% in the nuclear-profile function and l/m e in the electron- 
profile function. For pion production in the PrimEx experiment these cut-offs are numerically 
roughly the same, as follows from Eq. ( 1201) . 

At this point we draw special attention to Eq.()271) where we extracted the factor 
ki/k 2 \\ — E1/E2. Therefore, the Primakoff denominators, already encountered in the cross- 
section distribution of Eq.( fl6l) . take on a simple form when expressed in terms of transverse 
momenta, 

^l.J_ = ^ L_ (32 ) 

E 2 q 2 Q 2 k 2 2± + m 2 Q 2 ± + Qf 1 ' 

Remember that Q» is fixed, and does not depend on the transverse momentum variables. 
In Born approximation there are two amplitudes to be reckoned with; 

G(q,Q) = [F e (q)xF A (Q)].^, (33) 

#(q,Q) = F e (q)-F A (Q), (34) 

in terms of which the two-photon-fusion matrix becomes 

M 2l = iM 2l fc(q, Q) - #(q, Q) ia ■ h] \2M A y/E 1 E 2 ] , (35) 
e 4 

M 2l = Z g wll . (36) 

The expressions for the amplitudes F e (q) and F A (Q) are given in Eqs (|28|) and ( 123]) . In 
momentum space the nuclear amplitude is a direct product of amplitudes. In coordinate 
space, it is a convolution of potentials, 

G(q, Q) = J d 3 r e J d 3 r n e'^^^e"^^ [V e (r e - r ff ) x V A (iv)] • k u (37) 

with r e and iv the coordinates of electron and pion relative to the point nucleus. The 
argument of the electron potential is the relative coordinate between electron and pion. The 
definition of the spin-dependent function H(c\, Q) is similarly 

if (q, Q) = J d 3 r e J d 3 rv e^-^e-^ [V e (r e - r w ) • V A {r w )} . (38) 
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IV. COULOMB PHASE 



There are several improvements to the Born approximation that must be made. Among 
them we find effects of the elastic-Coulomb scattering of the electron, the distortion of the 
pion wave within the nucleus, and the extension of the nuclear-charge distribution. In this 
section we shall treat the elastic-Coulomb scattering of the electron. 

Both initial- and final-state electrons are relativistic and move along the eikonal in- 
direction. The Coulomb distortion of the electron wave function results in an eikonal-phase 



factor {3 1 

pixcQ>) — I 

v b r 

with b the relative impact parameter between the electron and the point nucleus, with a the 
cut-off radius of the Coulomb potential, and with 

T] = 2Za/v. (40) 

The velocity v can be put to unity, and the cut-off radius leads to an a-dependent phase 
factor common to all amplitudes. 

Changing the argument of the electron-profile function, replacing the electron-impact 
parameter b e by b e + b^, we get the Coulomb corrected amplitude 



e *oM = ('^T (39) 



G(q, Q) = J d%e^ J d%e^ 



_2o_T 
I b e + b n \_ 

x [r e (b e ) x T A (K)] ■ h- (41) 



The phase factor in the integrand of Eq. (1411) couples the integrations over the impact 
parameters b e and b^, and gives rise to a nucleus form factor, with a structure that could be 
quite complicated. It will not only depend on the absolute values of the momenta and q^, 
but also on their scalar product. Hence, it is expected to be significantly more complicated 
than the Coulomb-nucleus form factor encountered in pion- nucleus bremsstrahlung 6| . The 
phase variation must be properly investigated since Coulomb production through two-photon 
fusion interferes with coherent pion-nucleus photoproduction. The Coulomb phase of the 
coherent pion-photoproduction amplitude depends only on q± and is known analytically for 
point-like-nuclear charge distributions 

Unfortunately, in Eq. (T4TT) only one angular integration can be performed analytically. 
The three remaining ones must be done numerically. We start from the basic expression, 



Eq.flU 



G(q, Q) = ^ J d\e^ J dXe^ (b e x b w ) • & 



(2*) a 

■K 1 (m e b e )K 1 (Q\\b 1T 



2a 



hi 



(42) 



| b e + Wl 

Keeping the angular difference e — 0^ fixed, and integrating over the other angle gives as 
result 

m e Q\\ 

G(q,Q) 



, , . 6 e d6 e K x {m e b e ) \ bAK K^Q^) 
K^) Jo Jo 

/■27T 

•2tt y d</>[- sin(0 + 9 - <p Q )} ■ J ^ a/ (q±b e ) 2 + (Q±b w ) 2 - 2q ± b e Q ± b T cos i 



2a 



^/62 + bl + 2b e b n cos(0 + </> ? - 
The corresponding expression for H (q, Q) is obtained making the replacement 



sinfrA + 



— > COS(0 + 



(43) 



(44) 



Remember that in G(q, Q) and H (q, Q) the longitudinal component of q is defined to be 
q\\ = m e . 

We can check Eq. ([43l) by working out the Born approximation, i.e. the limit of vanish- 
ing Coulomb phase, rj = 0. Rewrite the sine function using the addition formula! Then, 
integration of the term proportional to sin0 obviously gives a vanishing contribution. For 
evaluation of the cosine term we apply the identity (see Appendix A) 



2tt 



2tt 



(45) 



d0cos(0)J o ^ V ' x 2 + V 2 — 2xy cos = 2-kJ 1 (x)2tt J\(y), 
to get as Born approximation 

fOO poo 

Cs(q, Q) = - sin(0 g - <f) Q )m e Q\\ / b e db e Ji(q±b e )K 1 (m e b e ) / fo^dfe^ Ji(Q±b v )K 1 (Q\\h K 

Jo Jo 



(q± x Q ± ) • fe a 







Q± 






[Q 2 \ 



(46) 



which, as expected, agrees with Eq. (!33|) . 

Before departing on a full-blown evaluation of Eq. (j4ip we take a look at the domain of 
integration. Since the modified Bessel function is exponentially damped for large values of 
its argument, there is effectively a cut-off in b e at the inverse electron mass; 



b ec = l/m e = 390 fm. 
10 



(47) 



The cut-off in b^ is at the inverse longitudinal momentum transfer Q\\. Its value depends on 
the meson produced. For the PrimEx experiment at 11 GeV/c, 



b vc { V ) = l/g||(77) = 14fm; 
6 TC (7r) = 1/Q||(tt) = 230 fm. 



(48) 



In the case of eta production the 5 e -region, with cut-off from Eq. (l47p . extends much further 
out from the nucleus than the 6^-region. As a consequence, it may be a good approximation 
to replace | b e + | by | b e | , thereby factorizing the integrand of Eq. (14 ip . All integrals can 
then be evaluated analytically, with the result 

Q±- 



G„(q, Q) = (q± x Q ± ) • h ^(ag) l V CT /i c (q 



Q 2 



(49) 



The definition of the point-Coulomb-form factor hc{c[), with q\\ = m e , is given in Eq.(21) of 
Ref.jfjj], where details of the integration can be found. In the same way the expression for 
the spin-flip amplitude becomes, 



ff„(q,Q) = qj_-Q± ^(ag) l V^ c (q 



Q± 



Q 2 



(50) 



f/,,(q, Q) differs from G v (q, Q) only in the extracted kinematic factor. There is in this 
approximation, effectively, only one Coulomb amplitude, and its g-dependent factor is also 
the Coulomb-form factor for meson photoproduction. 

We expect Eqs ( 1491) and ( 1501) to be useful for eta production, but for pion production 
there will be large corrections. But, since the two-photon-fusion amplitude .M27 interferes 
coherently with the photoproduction amplitude M.^ it is imperative to know the relative 
Coulomb phases of these amplitudes. In the photoproduction case, the virtual photon of 
four-momentum q is attached to a nucleon in the nucleus, and the radius of the nucleus is 
small compared with the cut-off b ec of Eq.( ]48~l) . The photoproduction amplitude will hence 
approximately factorize, one of the factors being the Coulomb scattering part of Eqs (1491 
and Q5D). 

Returning now to Eq. (j43|) . we observe that 



- sin(0 + d> a - 



-q_i_ • Qj_ sin + (q_i_ x Qj_) • k\ cos < 



(51) 



Integrating Eq. ( 143 j) with the sin factor gives a small return due to the approximate asym- 
metry of the integrand during inversion (f) — > —(f). The cos <j) term, on the contrary, gives a 
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large return due to the approximate symmetry of the integrand during the same inversion. 
We conclude that G(q, Q) is most important when <fi q — <pQ = ±7r/2, and vanishes identically 
when (j) q — (pg = 0, it, in concordance with the properties of G v (q, Q) of Eq.( 
Next, we observe that 



COS(0 + 



q_L ■ cos <p + (qj_ x Qj_) ■ k± sin < 



(52) 



By an argument similar to the one of the previous paragraph we claim that H (q, Q) is 
maximal when (p q — <J)q = 0, it and minimal, but not vanishing, when (p q — <J)q = ±7r/2. 

The two functions G(q, Q) and H (q, Q) are independent, but there are other ways of 
introducing independent form factors. We may, e.g., perform the 0- integrations of Eqs (143]) 
and (jSj) with factors of cos <fi and — sin 0. In Appendix B we show that precisely these 
integrands appear naturally if we remove the vector product (scalar for H) in G(q, Q) of 
Eq. fj4"Tl) . with the derivative trick. The new functions that emerge are labelled K(q, Q) and 
L(q, Q), and in terms of them 



G(q, Q) = K(q, Q) (q ± x Q ± ) ■ k x + L(q, Q) q ± ■ Q ± , 
H(q, Q) = K{q, Q) q ± • Q ± - L(q, Q) (q ± x Q ± ) • k v 

The integral definitions of the functions K(q, Q) and L(q, Q) are 



(53) 
(54) 



, . b e db e K^rriebe) / Kdb^ Kx{Q\\K) 
(2vr) 2 J Jo 



•2tt 



2- 



£(q, Q) 



•{COS0} 



m e Q 



■ Jo \ \/{.qJ>e) 2 + (Q±b 7T ) 2 - 2q ± b e Q ± b v cos< 
2a 

^b 2 e + bl + 2b e K cos(0 + <p q - 



11) 



(55) 



roo poo 

/ b e db e Ki(m e b e ) / b n db n K^Q^) 
(2vr) 2 Jo Jo 

2tt^ d0 • Jo (VfaiA) 2 + (Q±K) 2 - 2q ± b e Q ± b 7T cos<j?) 
{— sin0} 



2a 



v/6 2 + ft 2 + 2b e b n cos(0 + 



(56) 



In the Born approximation, r\ — 0, the L-function vanishes, L#(q, Q) = 0, whereas the 
i^-function is reduced to 

(57) 



**(q,Q) = ^- Q2- 
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a) |K/K/ b) |K/K B | 2 




123456 123456 
Q Q 



FIG. 2: Contour plots of the squared amplitude ratios for n production on Pb. In panels a, 
c, and d the ratios are with respect to Kq = l/(4Quqn), and in panel b with respect to the Born 
approximation A'^(q, Q). In panels a, b, and d the relative </>-angle is <p q — c/)q = vr/2, whereas in 
panel c it is <j) q — <j)Q = 0. The y-coordinate measures the momentum component q± in units of qu, 
and the x-coordinate the momentum component Q± in units of Qu. 

The numerical presentation of the functions K(q, Q) and L(q, Q) is difficult since they 
depend on three variables; q±, Q±, and the angle between them <p q — (pQ. We are attempting 
this in FigfSJ In Figs [2^, and [2j3 we present contour plots of the squared amplitude ratios for 
angle 4> q — 4>q = n/2. InEk the ratio is taken with respect to 

K = l/(AQ m ), (58) 

which is the maximal value of the Born amplitude fl57|) . attained at the double peak, and in 
|5b with respect to the Born amplitude itself. We conclude that in this configuration the K- 
amplitude vanishes, as does the Born approximation, when or q± approaches zero. Also, 
the magnitude of the i^-amplitude is considerably smaller than that of the Born amplitude, 
in the double-peak region, but approaches the magnitude of the Born amplitude when the 
values of Q± or q± becomes substantially larger than at the peak. 
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FIG. 3: The ratio -R(q, Q) defined in Eg . ([59]) graphed for pion and eta production. Note the dif- 
ference in scales! Each curve represents a fixed Q-value. The chosen values of the ratio Qx/Q\\ are, 
1.0 (solid line), 2.0 (dotted), and 4.0 (dashed). The g-values cover the range q±/m e = [0.1,10.0], 
with the largest g-value furthest away from unity. 

In FigfSb the relative angle is <p q — <pQ = 0, and we notice a radical change. The effect 
of the elastic Coulomb scattering of the electron is to make the i^-amplitude non- vanishing 
when Q± or q± vanishes. In fact, it has a maximum at Q± = q± = which is stronger than 
the Born maximum at the double peak. But when we go to larger values of Q± or q± the 
amplitude values are similar to those at <p q — <f)Q = ti/2. 

The L-amplitude vanishes when <f) q — <Pq = 0. Its relative magnitude for <ft q — (pg = n/2 
is graphed in FigfSJi. We again see the effect of elastic Coulomb scattering. The amplitude 
assumes its largest values when Q± and q± both take on values smaller than their peak 
values. But beyond the double peak the L-values are much smaller than the corresponding 
-fT-values, also for lead. 

Next, we return to the factorization property suggested in Eq. (j4"9l) . The function L(q, Q) 
is very small for eta production, since the angular dependence of the phase factor in Eq. ()55p 
may be ignored, leading to a vanishing integral. Hence, we are only interested in the K- 
function ratio 



i?(q,Q) = K(q,Q)/ ^(a#V CT Mq) • 



q- 



with the Coulomb form factor h c (q) defined as jg] 



h c (q) = T(2 - i77/2)r(l + irj/2)F(i V /2, 1 - i V /2; 2; ? ., ,. 

<Z_l + 9|| 



(59) 



(60) 



and q\\ = m e . Furthermore, we are only interested in situations where the variable Q± is in 
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the peak region. In Figj3] we have plotted the ratio -R(q, Q) for pion production in panel 
3a and for eta production in panel 3b. The curves traced correspond to fixed values of the 
ratio Q±/Q\\. The values chosen are 1.0, 2.0, and 4.0. Along the curves the q- values cover 
the range q±/m e = [0.1,10.0], with the smallest g-value, q = 0.1m e , associated with the 
points nearest to unity. We conclude that factorization is excellent for eta production. For 
pion production factorization is a good approximation for values of q in the peak region but 
becomes increasingly poor as the g-values increase. 
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V. EXTENDED NUCLEAR-CHARGE DISTRIBUTION 



The nuclear-electric field of a point-like-nuclear-charge distribution is given in Eq.(J2! 
and controls the photon radiation by the the nucleus. In the case of an extended nuclear- 
charge distribution Eq.( l2~4"j) is replaced by 

Vx(r) = ^ / dV|i^3 p(x'), (61) 

with charge distribution p(x) normalized to unity If we for simplicity consider a uniform- 
charge distribution with radius R u then Eq. (T2~4l) is replaced by 

-Is I 1, r > R u , 
V A (r) = —^{ (62) 



The profile function is implicitely defined in Eq. (12511 and as a consequence 

o poo 

T A (b, Q||) = — / dz e- iQ * z V A (b, z). (63) 

When b > R u the trajectory always runs outside the nucleus, and the expression for the 
point-like case applies, Eqs (12~6|) and (l62l . 



b > Ru : r A (b,Qn) = ^ [6g,|ifi(6Q||)]. (64) 
The case b < R u is more complicated and integrals cannot be done analytically; 

b<R u : r A (b,g,|) = \ [-^^^ si n(Q|| V ^7^) + 6Q||^ 1 (6Q||) 



qb 



Q\\Ru R^ 



-by/l-P/E*-b ^ dz (62+ "^ 2)3/2 j- (65) 

Here, the first term inside the brackets is the contribution from the z-integration over the 
uniform interior. The second term is the total contribution for a point-like charge distribu- 
tion, and the last two terms removes from it the contribution to the ^-integration over the 
nuclear interior with a point-like charge distribution. 

In the production of eta mesons we expect a substantial part to come from the nuclear 
interior, and hence, a proper description of the nuclear-charge distribution is certainly called 
for. For production of pi mesons, on the other hand, very little comes from the nuclear 
interior and it should be sufficient to treat the nuclear-charge distribution as point like. 
Since, we are mostly interested in the pion case we do not present any numerical results for 
extended-charge distributions. 

The profile function of the electron is not affected by this discussion. 
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VI. PION MULTIPLE SCATTERING 



The pion is produced at r = (b, z) and moves from there with constant b to z = oo. 
Along this trajectory it can scatter off the nucleons in the nucleus. Following Glauber [s| 
this scattering introduces a pion- wave-function distortion represented by the factor P(r), 
which for heavy nuclei takes the form 

P(r) = exp[-^T(b,z)], (66) 
with the target-thickness function 

POO 

T(b,z) = / dz'n(b,z'), (67) 

J z 

and n(r) the nucleon density in the nucleus. This density is normalized to A, the number of 
nucleons in the nucleus. For the structure of the density, we refer to Eq. (ll25j) of Appendix 
C. 

The integration along the path of the pion gives for the uniform nucleon-density distri- 
bution, 

'0; b > R u 

T{b,z) = { , (68) 

n^R\ -b 2 -z; b < R u , \z\ < ^ R 2 U - V 

0; b<R u , z> y/B* - 6 2 

with n the nucleon density. 

The pion rescattering is taken into account by replacing the nucleus potential of Eq. §M§ 

by 

v " (r) = 4^FT^F exp[ -y T(b ' 2)] - (69) 

with r = (h,z). This extension could have some influence on eta production, since a large 
fraction of this production takes place inside or near the nucleus. For pion production it 
should be negligible. 
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VII. COHERENT-PION PHOTOPRODUCTION 



The neutral pion can also be produced in a coherent-nuclear-photoproduction process, a 
contribution which interferes coherently with the two-gamma-fusion process. Since we are 
interested in a coherent-nuclear-reaction amplitude its main contribution is generated by 
the spin-isospin non-flip term of the underlying nucleon amplitude. This nucleon amplitude 
can be modelled by the w-exchange diagram of Fig. HI In writing down the contribution of 
this diagram we are particularly anxious to get the correct Lorentz structure. 




FIG. 4: Graph describing pion-photoproduction on the nucleon through i-cannel omega exchange. 
The cu-decay vertex into 7r°7 in the graph of Fig. H] is described by the matrix element 



(70) 



where we may replace by p n . The decay width of T(cj — > 7r°7) = 760 keV sets the coupling 



strength to g { 



0.322. 



The vetex for uj coupling to a nucleon takes the form 



M(uNN) = -ig U NNu{p') fa + — %- % tJ) u(p). 



(71) 



Numerical values for the coupling constants are g U NN — 15.9 and = 0. Furthermore, in 
coherent-nuclear production the nucleon spin-flip amplitudes are heavily supressed. There- 
fore, if we rewrite Eq. (|7ip with the aid of the Gordon decomposition 



u(p), 



(72) 



it follows that only the current part survives in the nuclear matrix element. The effec- 
tive nucleon-photoproduction amplitude responsible for the corresponding coherent-nuclear 
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amplitude becomes 

M{lN n o N) = ^^1^)^^^+^. (73) 

The matrix element for ir° photoproduction is equal for protons and neutrons. 

In order to make this expression useful at higher energies we replace the omega-pole factor 

P u = r-^ (74) 
t — m£j 

by its Reggeized version i.e., 

RM,t)={-) . ™» .-J—.±L±±^ . (75) 

\8 J sm(7ra w (t)) r(a w (t)) 2 

The signature of the omega trajectory is positive, = 1, and its parametrization has been 

determined from photoproduction experiments by Guidal et al. Q to be 

a u (t) = 0.44 + 0.9*, (76) 

with * in units of (GeV/c) 2 and so that a'^ = 0.9 (GeV/c) 2 . In this study a strong emphasis 
is on reproducing the minimum in the cross-section distribution at — * = 0.5 (GeV/c) 2 , 
whereas in our nuclear application the relevant t-values are substantially smaller. 

Specializing to the lab. system where the initial proton is at rest, we get as effective 
nucleon-photoproduction matrix 

M L (jN -> 7r°iV) = ie^g wNN R w (s, t) (k 7 x e 7 ) ■ K [2M] . (77) 

As we can see there is no spin-dependent contribution. For insertion into the coherent- 
nuclear-photoproduction amplitude only the forward- angle-pion photoproduction amplitude 
is needed, so we can set * = 0. Furthermore, we have chosen to present numerical results 
for photon energies of k n = 10.7 GeV/c, corresponding to s = 21 (GeV/c) 2 . At this energy, 
and at * = (GeV/c) 2 , 

L u (s, t) = ^g^NNR^s, t) = 5.86 - *4.85 (GeV/c)" 3 . (78) 

Also rho exchange contributes to neutral-pion photoproduction, but with amplitudes of 
opposite signs for proton and neutron. Taking the parameters of the rho trajectory again 
from Ref.|7J, and a rho-decay width of T(p — > 7r°7) = 68 keV, we get at s — 21 (GeV/c) 2 
and * = (GeV/c) 2 , 



m-n 



L p p (s, t) = -L n p (s, t) = ^g pNN R p ( Sj t) = 0.169 - z0.198 (GeV/c)- 3 . (79) 
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The rho parameters are much smaller than the omega parameters, in part because the 
coupling constant g P NN = 3.25 is much smaller than the corresponding omega-coupling 
constant. In addition when we calculate the nuclear amplitude the omega contributions 
from protons and neutrons add, yielding a multiplicative factor Z + N = A, whereas the 
rho contributions subtract, yielding a smaller multiplicative factor Z — N. Consequently, it 
is legitimate to ignore the rho contribution altogether, and we shall do so. 



fa k 2 




Pi P2 

FIG. 5: Graph describing coherent-pion photoproduction in the Coulomb field of a nucleus in 
inelastic-electron coherent-nucleus scattering. 

The nuclear-photoproduction diagram is shown in Fig. [51 The external kinematics is the 
same as in Fig. (TJ and the electron vertex the same as in Eq. (JTj) . The photon propagator is 
also the same, but for the photoproduction vertex we need the nuclear version of Eq. fl73l) . 
However, the similarity with two-photon fusion does not end here. The coupling of the 
omega to the nucleon, Eq.( 17TT) . is the same as the coupling of the photon to the proton, 
Eq.(t3), provided we neglect the magnetic terms, which we do. In Born approximation the 
photoproduction amplitude has thus the same structure as the two-photon-fusion amplitude, 
Eq. ffT3l) . Before Reggeization, the photoproduction amplitude is obtained by changing the 
photon-nucleus coupling constant Ze into the omega- nucleus coupling constant Ag^NN, and 
the photon propagator 1/Q 2 into the omega propagator 1/{Q 2 — m^)- 

But there is also a nuclear form factor. From Eq. (1191) we know that the momentum 
transfer to the nucleus along the direction k\ of the incident electron is fixed 

2 

a 77? 

Qll = (K - q) • h = (80) 

whereas the component orthogonal to this direction varies. The nuclear form factor 
must also reflect the distortion of the pion wave due to hadronic scattering within the 
nucleus. Moreover, the nucleon amplitude of Eq.( 177|) vanishes in the forward direction, as it 
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proportional to (q x 



Q. The mathematical form of the appropriate nuclear form factor 



F S (Q) is given in Ref.|5|, and its operational definition is through 

/poo 
d 3 xe- 4Q - x exp[-i<y &z'n{h,z')\ h- V 6 n(b,z), (81) 

where h = q x e 7 is a vector in the impact plane. The pion is produced at the point (b, z) in 
the nucleus where the nucleon density is n(b, z). The incoming photon wave is undistorted 
whereas the distortion of the outgoing pion wave is determined by the pion-nucleon total 
cross section a'^ = a n (l — ia w ). 

To explain this result we start from the omega-exchange factor 

h ' Q / d 3 re" lQ r (-ih ■ V^—e~ m - r j . (82) 



Q 2 + m 2 J \ Airr 

Then we fold with the nucleon density n(r) and take the limit of large omega mass, 

J d 3 re" iQr J d 3 r N (-ih • ^ T ^-^— e -™»\*-*s\\ n ^ N ) 

= [ d 3 re- iQ - r [ d 3 r N —^ re"^ 1 ™ 1 (-zh ■ V rjv n(rjv)) 

J J 47r|r-rAr| 

= At / d 3 re- lQ ' r (-ih • V 6 n(r)) . (83) 
™>i J 

Adding distortion due to pion scattering leads to Eq. flgTj) . In this derivation we have used 
the standard expression for the omega pole instead of its Reggeized version. This should 
not change too much since we essentially only need the forward value of the nuclear form 
factor. 

All this adds up to a nuclear photoproduction matrix element, which in Born approxi- 
mation is given by 

MJe-A -)• e-7i A) = M^(e~ -»• 1 e~)^^M v hA -> ir°A) 

= -i(ie)(ze)L w (<)F s (Q)^ u{k 2 )Yu{k x ) e^ p ^K"{p l + p 2 ) CT ,(84) 

where s 7 is the energy of the underlying nucleon photoproduction process. In our application 
s 7 = 21 (GeV/c) 2 . Evaluating the matrix element in the lab. system gives a formula similar 
to the one for the two-gamma matrix element, Eq. (1131) . 

M y = ie 2 2M A L u { Sj )F s {Q) 1 u(k 2 )j • (q x Q)u{h), (85) 
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where we have replaced K by Q. 

We have previously shown that in the impact parameter plane the integration over the 
radiation potential extends out to 396 fm. The extension in impact parameter space of the 
photoproduction amplitude is limited to the nuclear region, which never extends more than 
7 fm from the origin. In the folding of the two potentials we may neglect the extension of 
the nucleus and factorize the two amplitudes. Thus, to a good approximation the electron 
multiple scattering is contained in the expression 

M-f = iM 7 M(q, Q) (q_L x Q ± ) • h - q ± • Q ± icr ■ h 2M A ^/E 1 E 2 \ , (86) 

where ftf 7 represents the product of coupling constants 

A/; = e 2 iazig^^s, 0KF s (Q), (87) 

and 

M{q, Q ) = — j 5 h c (q). 88 

ml q 2 

This amplitude includes multiple-pion scattering through the nuclear form factor F S (Q). In 
the present application we may evaluate F S (Q) at Q = 0. For a uniform nuclear-density 
distribution an analytic expression for F s (0) can be obtained. This is demonstrated in 
Appendix C. Also, we should keep in mind that in the form factor hc{q) the argument is 
interpreted as q = (q_i_,m e ). 

We have defined the function M(q, Q) of Eq.( l88|) so that its origin in omega exchange 
becomes obvious, since ml—Q 2 ~ m^. The dimensionless factor Rcj(s, 0)m^ then determines 
the strength of the omega-exchange pole as obtained from the Regge-pole model, and 

0)ml = 0.095 - i0.078. (89) 
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VIII. CROSS-SECTION DISTRIBUTIONS 



The cross section distribution is given by Eq.(|2]) of Sect.l, 
da 1 



d 2 k 2± d 2 K ± dk n 32(2n) 5 k 1 E(k 2 )K ll Ml 



\M 2l + M 



71 ' 



(90) 



with K\\ the component of pion momentum along the incident direction. The expression 
for the two-photon fusion matrix is found in Eqs (1351) and ( )36l) 



M 2l 

M 2l 



2~ 



G(q, Q) - #(q, Q) to ■ fcj \2M A ^/E^E : 



(91) 
(92) 



The functions G(q, Q) and i?(q, Q) are defined in Eqs (1331) and flSJ), but we shall rather 
employ their decompositions into functions K(q, Q) and L(q, Q) of Eqs ( )53l) and ( 1541) . 
The expression for the photoproduction matrix is found in Eqs fl86|) and (jSj 



M 7 



M 1 M(q, Q) [(q ± x Q±) ■ h - q± ■ Q± icr ■ k] ^NlA^/EhE, 



2 doJify 



gwNNRu>(s, 0)m 2 J F s (0), 



(93) 
(94) 



with the function M(q, Q) as defined in Eq. fl88|) . In the peak region the form factor F S (Q,) 
can be evaluated at zero momentum transfer, where it gives the effective number of nucleons 
in omega exchange. According to Appendix C we have for lead, with pion rescattering 
included, F a {0) = 0.70 A. 

From this exposition we derive for the unpolarized cross section the expression 



dcx 



d 2 k 2± d 2 K ± dk 2 



7rK\ 



K(q,Q)-A^M(q,Q) + L(q, Q) 



2 n 



(95) 



with Mr the ratio of coupling constants of the photoproduction and the two-photon fusion 
processes, 



Mi 



R 



9unT-y9uiNN 



RumlFsW/Z. 



(96) 



The expression for Mr is obvious, considering the similarities of the two processes, provided 
we remember that -R^m 2 is the strength of the omega pole, according to the Regge modell. 
An interesting feature of Eq.( l95l) is that the photoproduction M-amplitude only interferes 
with the two-photon-fusion ^-amplitude. 
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First, we would like to determine the relative size of the photoproduction contribution. 
We do this in the Born approximation where 



**(q,Q) = S'^> (97) 



q± Qj 
q 2 Q 2 

L B (q,Q) = 0, (98) 

M B (q,Q) = %.%. (99) 
q 2 ml 

In this approximation the g-dependence is the same for the two-photon-fusion and the pho- 
toproduction amplitudes. The ratio between the coupling-constants factors is 

\J\f R \ = 3.08 x 10" 3 . (100) 

For large values of Q the photoproduction contribution dominates. Only in the vicinity 
of the peak can the two-photon fusion contribution be measured. At the peak, where 
Q± — Q\\ — 0-85 MeV/c, the enhancement factor becomes 

4.22 x10 s , (101) 

making the two-photon-fusion amplitude 1300 times stronger than the photoproduction 
amplitude. Consequently, at the peak a clean measurement of the pion-decay constant 
should be possible. This conclusion is not affected by our ignorance of the relative sign 
between J\f 2y and N^. 

Even if the two-photon-fusion contribution is dominant in the vicinity of the peak, the 
theoretical analysis is non-trivial. There are two contributions represented by the K- and 
//-amplitudes. They both depend on the relative angle between and q^. Only the K- 
amplitude interferes with the photoproduction amplitude. The i^-amplitude is dominant 
at all angles but is considerably smaller than its Born approximation. The L-amplitude 
vanishes in the Born approximation, but is for lead at some angles as important as the K- 
amplitude and at other angles unimportant. So the theoretical description of pion production 
in the peak region is certainly non-trivial. A closer scrutiny of Fig. 2 emphasizes this picture. 

If we aim for a determination of the pion decay constant g nll , it is mandatory that the 
transverse momentum component of the pion, Q±, be in the peak region, Q± ps Q\\, but 
the transverse momentum of the electron, q±, need not be in its peak region, q± m e . 
The point, however, is that if we are in this peak region as well the cross section is at its 
strongest. 
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IX. APPENDIX A 



The integral encountered in Eq.( l45|) is a special case of an integral by Sonine and Gegen- 
bauer jsj. It is most easily established by expanding the Bessel function using Neumann's 
(or Gegenbauer's) addition theorem. Since the integral is not always listed in standard 



tables of integrals, e.g. missing in Gradshteyn and Ryzhik Jq], we shall here show how to 
evaluate it. The method is simple and employs the same technique that is used for other 
impact-parameter integrations of Sect. IV. 
Define the function F(q, Q) as 

F(q,Q)=/ dtp a d(p b (a x b) • hexp[iq- a-iQ ■ b]. (102) 
Jo Jo 

Here, the vectors q, Q, a, and b, are all in the xy- impact plane and h along the positive 
z-direction, i.e. along the normal to the impact plane. Integrations are over the angles p> a 
and (fib. The steps of integration are as follows 

I /*27r /*27r 

■F(q> Q) = — (V 9 x V Q )-h J dcp a J d(p b exp[iq- a - zQ • b] 

= ^(V 9 x V Q ) • h 27rJ {qa)27rJ {Qb) 
= 2irJ 1 (qa)2'KJ 1 {Qb) (qxQ)-h 

= sin(<y9g — ifq) 2ir.Ji(qa)2nJi(Qb). (103) 

The second more complicated approach is that of Sect. 3. We keep the angle ip = ip a — (p h 
fixed while integrating over (p b . During the integration over (p^, the factor 

(a x b) • h = -sinv? (104) 

stays constant as well as the scalar product a-b = cos ip. The exponential can be manipulated 
into 

E = q- a- Q ■ b = VAcos(ip A + <Pb), (105) 

with 

A = (qaf + (Qb) 2 - 2qaQbcos(<Pg - <Pq - (p), (106) 
and the angle ipA defined by 

, qasinjipg - <p) - Qbsm{<p Q ) 

tan(^ A = r-s r — — — -, (107) 

qacos(<p q - <p) - Qbcos{ip Q ) 
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a knowledge which is not really needed. After integration over we get 

I-2tT 

F(q, Q) = 2tt / d<p[- sm{cp + <p q - <p Q )]J (y/A(vj), (108) 
Jo 

where now 

A(<p) = (qa) 2 + (Qb) 2 - 2qaQb cos (p. (109) 
Furthermore, from the angle-addition theorem 

sm(ip + ipg- v?q) = sin(^) cos(y? 9 - y?g) + cos(^) sin(y2 g - <^q), (110) 

we may conclude that the term proportional to simp vanishes on integration. Comparing 
with expression (11031) it follows that 

2tt J d0cos(0)J o (\A 2 + y 2 - 2xy cos <f?J = 2vr J 1 (x)2n J x (y) , (111) 

as promised. If we need integrals of this type but with higher powers of cos <fi in the integrand 
we start from Eq. (jl02j) with integrands containing higher powers of (a x b) ■ h. 
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X. APPENDIX B 



The alternative approach to the nuclear form factor of Eq. (HIT) is to replace b e and by 
derivatives 



b e -»■ -iV, 



— z 



q— + (k x g)-^— 



(112) 
(113) 



Here, q and Q are two-dimensional vectors in the impact parameter plane, and k is along the 
normal to this plane, i.e. along the positive ^-direction. As a consequence the polar-angular 
basis vectors are 



'H>Q 



k x q 
k x Q. 



In the evaluation of Eq.( HT|) and its counterpart we meet the differential operators 



V 



A 



d d 1 d d 
+ 



V, 



dq dQ qQ dcp q d(p Q ' 
1 d d 1 d d 



q dip q dQ Q dip Q dq ' 



in terms of which 



b e ■ b n V, ■ V Q = q ■ Q V K + k- (q x Q) V L , 



k-(b e xb n ) ->• • (Vg x Vq) — k - (q x Q) V k — q ■ Q V l . 



(114) 
(115) 

(116) 
(117) 

(118) 
(119) 



The basic functions G(q, Q) and H(q, Q) are replaced by two other functions K(q, Q) 
and L(q, Q), where L(q, Q) vanishes when the elastic Coulomb is turned off. 
Going back to Eq. (1431) . once more, we define the new form factor if (q, Q) as 



K(q,Q) 



j m e Q 



/•27T 

2tt y d0- Jo (V^e) 2 + {QKY - 2qb e QK cos. 



2a 



- vq 



^ b 2 + b 2 + 2^ cos(0 + 



(120) 



For the form factor L(q, Q) we substitute the operator T>l for T>k- Finally, remember that 
in this appendix q and Q represent the transverse momentum components. 
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Pulling out a differential operator from inside an integral is a useful approach when the 
remaining integral can be done by hand. That is not the case here but the technique nev- 
ertheless gives new insights. We start with L(q, Q). In differentiating the elastic-Coulomb- 
phase factor with respect to an angle we note that d Vq = —d VQ = d v . Then we make a 
partial integration and let act on the Bessel function. The operator T>l now acts on the 
Bessel function, and with 

X = (qb e f + {QKf - 2qb e Qb n cos <f) (121) 

we get 

V L J (VX) = ^d^qd q + Qd Q )J (VX) = J (VX). (122) 

The differentiation of the Bessel function involves applying the operator qd q + Q8q to a 
homogeneous second order polynomial, X, an operation resulting in a factor of two. 
Repeating these operations for K(q, Q) we get 

VkMVX) = (d q d Q - dA Jo(VX) = -6Acos0 J (VX). (123) 

After reintroducing our earlier notation for the transverse components, the integral defi- 
nitions of the functions K(q, Q) and L(q, Q) become as spelled out in Eqs (1551) and (1561) . 
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XI. APPENDIX C 



The pion-hadronic form factor F S (Q) of Eq.( )8T|) can be expressed as 



F S (Q) 



-2tt I dze 

-oo 



bab — - — : — exp 



1 / 



dz'n(b, z') 



(124) 



In the Coulomb region, the variation with momentum transfer of this form factor is weak. 
Therefore, it is sufficient to consider the forward direction, setting Q = 0. For a uniform 
nuclear density, the nucleon density is 



n(r) = n 9(Ru - r), 

with R u = tqA 1 ^ the uniform nuclear radius, and the central density 

3 3A 



(125) 



n 



If we observe that 



dn(r) 



-n -5(r - R u ), 



(126) 



(127) 



db r 

the integrations in Eq. (11241) could easily be done by hand. After a few straightforward 
manipulations we get 



with 



C = a' n n R u - 



128) 



(129) 



A typical value of the pion-nucleon total cross section is a n = 26 mb. For a n = the 
numerical value of the form factor becomes, for lead, F s (0) = 0.70 A. In the absence of pion 
scattering, ( = 0, and F s (0) = A. 
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